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Abstract 

We show how a central limit theorem for Poisson model random polygons implies a 
central limit theorem for uniform model random polygons. To prove this implication, 
it suffices to show that in the two models, the variables in question have asymptotically 
the same expectation and variance. We use integral geometric expressions for these 
expectations and variances to reduce the desired estimates to the convergence (1 + 
— ^ e" as n — ^ CO. 



1 Introduction 

Given a convex set C of unit area, we may define two random polygon models. For the 
Poisson model, we consider a Poisson process of intensity A inside K, and define ^ to be 
the convex hull of the points of this process. For the uniform model, we take n independent 
random points distributed uniformly in K and let PK,n be their convex hull. For a polygon 
V inside K, we let N{V) denote the number of vertices of V, and we let A{V) denote the 
area of K\V. When one wants to prove central limit theorems for and A in either of the 
two models of random polygons, it is often the case that the Poisson model is easier than 
the uniform model. Indeed, in general, results are proved first for the Poisson model, and 
then more arguments are needed to deduce a corresponding result for the uniform model. 

Recently, the author [1] studied the Poisson model of random polygons and proved the 
following central limit theorem for and A: 

Theorem 1.1 ([4J). As X oo, the following estimates for IIk,\ hold uniformly over all K 
of unit area: 



sup 
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Here = P{Z < x) where Z is the standard normal distribution. 

In this paper, our goal is to show how to derive the following corollary for the uniform 
model: 
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Corollary 1.2. As n ^ oo, the following estimates for Px.n hold uniformly over all K of 
unit area: 

sup P ^ ^ <x \ - $ X 



(1.3) 



sup 



^ (1.4) 



Here $(a;) = P{Z < x) where Z is the standard normal distribution. 

From the estimates derived in this paper, a secondary result from |3] (Theorem II. 3l below) 
also carries over immediately to the uniform model. We should say that Theorem 11.31 and 
the consequence derived here, Corollary 11.41 have both been proven independently by Imre 
Barany and Matthias Reitzner. 

Theorem 1.3 ([4J). As \ oo, the following estimates for TIk,\ hold uniformly over all K 
of unit area: 

E[A^] X Var X AE[A] x Var A (1.5) 

Corollary 1.4. As n oo, the following estimates for PK,n hold uniformly over all K of 
unit area: 

¥.[N] X Var AT X nE[A] x Var A (1.6) 

Though both Theorem 11.11 and Corollary 11.21 have been known for quite some time in 
the case that either is a polygon or dK is of class C^, the proof given here of Theorem 
11.11 =^ Corollary 11.21 for arbitrary convex K appears to be new. Corollary 11.21 answers a 
question of Van Vu [Ij . 

We will see below that in order to prove that Theorem 11.11 implies Corollary 1 1.21 it suffices 
to show that when n = A, the random variables iV(Px,n) and NiJlK,\) (as well as A{PK,n) 
and A{I{k,\)) have the same expectation and variance up to a small enough error. This 
is essentially the same strategy used by Van Vu in [7j to derive a similar implication for a 
special case of random polytopes. 

However, in contrast to [7j, we will use relatively down to earth integral geometry to 
establish our estimates, instead of sophisticated arguments from probability theory. The 
approach we take is conceptually very simple. We write down integral geometric expressions 
for the expectation and variance for the Poisson and uniform models, and then estimate 
their difference in the limit n = A — > oo. In this formulation, the "reason" that the desired 
convergence holds is completely transparent: it is essentially reduced to the convergence 
(1 + — )■ as n — )■ oo. Also, the variables A^ and A are treated simultaneous with an 
identical proof for each (c.f. where the case of /j, the number of i-simplices, is harder 
than the case of the volume and requires a new idea). An admitted disadvantage of this 
approach is that one has to actually write down these integrals explicitly, however once this is 
done, no further manipulations are necessary. It is interesting to observe that using integral 
geometry is almost never the "right" way to prove statements along the lines of Theorem ll.il 
or even Theorem II. 3[ essentially because the expressions quickly become too complicated to 
deal with either conceptually or theoretically. However for our applications here, the desired 
estimates become simple when written in terms of the integral geometry, so we in fact believe 
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that these integral geometric expressions do in some sense give the "right" proof of our main 
lemmas. 

One expects that our results and the proofs given here will admit straightfoward gener- 
alization to higher dimensions. For random polytopes in dimension d > 3, Theorems II . lf[L^ 
are all known in the case of fixed K whose boundary is and has nonvanishing Gauss 
curvature, due to Reitzner [5] and Vu [7]. In the case that is a polytope, the analogue 
of Theorem 11.11 was proven very recently by Barany and Reitzner [5] (one expects that an 
analogue of Theorem 11.31 also follows from their methods). We expect that if applied to 
higher dimensions, the methods in this paper would show that Corollaries 11.21 and 11.41 follow 
in any situation in which Theorems 11.11 and 11.31 respectively are known to hold (in particular 
for the case that if is a polytope). It is conjectured that Theorems ll.lflL^ hold for d > 3 
with no restriction on K. 

1.1 Acknowledgement 

We thank the referee for useful comments and in particular for making us realize an error in 
the original proof of Lemma 13.11 

2 Notation and definitions 

We now review some definitions and two basic lemmas from 

In this paper, K will always denote a (bounded) convex set in M^. Any constants implied 
by the symbols or x are absolute; in particular they are not allowed to depend on 

K. 




(a) Illustration of W{e) (b) Illustration of Ck{p, 0) 



Figure 2.1: Illustration of some definitions. 

Many of the following definitions are illustrated in Figure 12.11 We may leave out the 
subscript K later when doing so is unambiguous. 

Definition 2.1. We define the random variable W'p{9) to be the vertex of V which has an 
oriented tangent line at angle 9. This is illustrated in Figure [2?]T a). 

Definition 2.2. A cap at angle 9 is the intersection of K with a halfplane Hq at angle 9. 
We may specify a cap at angle 9 by giving either its area r or a point p G dHe. These are 
denoted CK{r,9) and Ck{p,9) respectively; the latter is illustrated in Figure l2T]T b). 
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Definition 2.3. We define the real number Ak{p, 0) to be the area of the cap Ck{p, 0). 

Lemma 2.4. The random variable Wii^x{9) has probability distribution given by A exp{—XAx{p, 0)) dp 
where dp is the Lebesgue measure. This has total mass 1 — e-^Area(x)^ ^ empty with 

probability e-^Area(_ft')_ 

Proof. This follows directly from the definition of a Poisson point process. The probability 
that no point lands in Ck{p,G) is exp{—XAK{p,0)), and we multiply this by Xdp, which is 
the density of the Poisson point process. 

Alternatively, we may differentiate exp^—XA^ip, 0)) with respect to the direction orthog- 
onal to 6 and divide by the length of OHq fl K. This also yields A exp(— Ay4x(p, 0)) dp. □ 

Definition 2.5. We define the function fxix^O) : [0, 1] x ]R/27r — > M as follows: 

( length of (dHe) n K where Ci^(log ^,e) = HgnK 
fK{x, 9) = I iix> exp(- Area(K)) (2.1) 

[o if X < exp(- Area(ii')) 

It will be important to have the following bound on the growth of /: 

Lemma 2.6. If y < x, then: 

fjy) < fix) 2) 
V-logy ~ V- log X 

The bound above is sharp, for instance /(x) = const ■^/^^Togx for K = {x,y > 0} (i.e. 
the first quadrant). 




angle 6 



Figure 2.2: Illustration of the function h. 



Proof. Project K along the lines at angle 9 to get a height function h : [0, oo) — )■ M>o; 
in Figure 12. 2[ h{£) is the length of the thick segment. Now if A{i) = h{£') di' then 
f{exp{—A{i))) = h{i). Thus we see that it suffices to show that the function: 

h{i) 



(2.3) 
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is decreasing. Differentiating with respect to we see that it suffices to show that: 

h{ey - 2h'{i)A{i) > (2.4) 

For i = 0, the left hand side is clearly nonnegative, and the derivative of the left hand side 
equals —2h"{i)A{i), which is > by concavity of h. □ 

When proving central limit theorems, it is important to decompose N and A into local 
pieces. Thus we define N{a, (3) to equal the number of edges with angle in the interval {a, (3). 
Then it is easy to see that: 

N = N{ai, as) + N{a2, 02) + ■ ■ ■ + N{aL, ai) (2.5) 
A similar decomposition is valid for A, where A{a, (5) is best explained graphically in Figure 




Figure 2.3: Illustration oi A{a,(3) 



Consider for the moment the Poisson model, and let X denote N 01 A. In [3], it is 
shown that if one chooses the partition so that each interval [a^, aj+i] has constant affine 
invariant measure (a notion from [1] which will not concern us here), then X{ai,ai^i) has 
constant expectation and variance, and the correlation between X{ai, Oj+i) and X{aj, dj+i) 
is exponentially decreasing in \i—j\ (specifically, an a-mixing estimate is proved). From these 
facts, along with a general lower bound on the variance of X due to Barany and Reitzner [2] 
(their result also holds in higher dimensions), it follows on general principles that a central 
limit theorem holds for X in the Poisson model. 

3 Proofs 

Let us begin by stating the lemmas which we will prove. 
Lemma 3.1. As n ^ 00, we have: 

sup \P{A{UK,n) <X)- P{A{PK,n) < x) | ^ (3.1) 
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sup \P{N{UK,n) <X)- P{N{PK,n) < x) | ^ (3.2) 

X 

uniformly over all convex K of unit area. 
Lemma 3.2. As n ^ oo, we have: 

mN{UK,n)]-nN{PK,n)]\ = o{^VaT N {UK,n)) (3.3) 
\E[A{UK,n)] - m{PK,n)] \ = o{^JVaT A{UK,n)) (3.4) 

uniformly over all convex K of unit area. 

Lemma 3.3. As n ^ oo, we have: 

Var iV(n^,„) ~ Var N{PK,n) (3.5) 

Var A{I{K,n) ~ Var A{PK,n) (3.6) 
uniformly over all convex K of unit area. 

Lemma [XT] says essentially that as n — )■ oo, the functionals of PK,n and Il.K,n have asymp- 
totically the same distributions. Lemmas 13.21 and 13.31 give us equivalence of the expectation 
and variance in the two models. It is elementary to observe that these lemmas combine to 
give Corollary IL2I (and also Corollary II. 4p . We note that Van Vu has observed (13. 4p in [TJ 
p224, Proposition 3.1 and Remark 3.5]; our proof is different. 

3.1 Proof of Lemma 13JJ 

Proof. We follow and slightly correct an argument of Reitzner [5l p492-3]. We thank the 
referee for asking us to clarify our use of Reitzner 's argument, since it was this that led 
us to realize the error. The argument below becomes fallacious if we write P{X < x\E) 
(as Reitzner does) everywhere we have P{X < xhE). The problem is that the the first 
equality in equation (13. 9 p (which is Reitzner's equation (13)) is false in this case (note our 
E is Reitzner's A). Thus the proof in [5] is typographically very close to being correct; there 
is no problem once we replace every P{X < x\E) with P{X < xhE). 

First, fix e > and let = UeGR/27r ^(^' ^) union of all caps of area e. Let Ep 

and Eyi be the events that dP and dU are completely contained in respectively. Now 
trivially, we have that P{Ep), P{Eji) — )■ 1 as n — )■ oo, uniformly over all K of unit area. 

If B is any event, then \P{B) - P{BkE)\ < 1 - P{E). Letting 5 be X < x, we have: 

sup sup \P{X{IiK,n) <x)- P{X{IiK,n) < X & En) | ^ (3.7) 
K X 

supsup|P(X(P;^,„) <x)-P(X(P^,„) <x&Ep)| ^0 (3.8) 

K X 

as n ^ OO, where X denotes either N or A. 

Now consider P(X(ni^,„) < xkE^) and P(X(Pk,„) < xkEp). Observe that if we 
condition both probabilities on the number of points of the process in 5"^, then they become 
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equal. Let us call this probability P{X < x Sz E\k). In other words, suppose we place k points 
uniformly at random in S^. Then P{X < x^E\k) is defined to equal the probability that 
the boundary of their convex hull is contained in and X < x. Thus setting p = Area(S'e), 
we find: 



\P{X(IIk,„) <xkEn)- P{X(Pk.„) <xkEp) 



E 

fc=0 



.^-npp(^X <xkE\k) 



n 



< 



E 

A:=0 



-np 



pY-''P{X < xkE\k) 

n—k 



k\ 



n 



p^[l 



pY 



< 2p (3.9) 



where (^) = if k > n. The last bound is due to Vervaat [B]. Combining fl3.7p and fl3.8p 
with (13. 9p . we find that: 



lim sup sup sup \P{X(JlK,n) < a;) — P{X{PK,n) ^ x)\ < 2sup Area(S'£ 

rt— >oo K X K 

But we may choose e > arbitrarily, so we are done. 



(3.10) 



□ 



It may indeed be possible to take a similar strategy to prove Lemmas 13.21 and 13.31 How- 
ever, in this case bounding the sum in equation (13. 9p becomes harder, since we have ex- 
pectations instead of probabilities, and the former are not bounded by 1. Also, proving 
analogues of equations (13. 7p and (13. 8p becomes nontrivial. Since we need good estimates 
for Lemmas 13.21 and 13. 3[ choosing e correctly as a function of n and estimating P{Eyi) and 
P{Ep) becomes an issue. 



3.2 Proofs of Lemmas [33] and [3731 

The proofs of Lemmas 13.21 and 13.31 will make use of some simple integral geometric expres- 
sions for the expectations and variances in question. The derivation of these expressions 
is completely elementary. The integrals appear complicated, though the point is not their 
exact form, but rather that they are almost identical for Pk and 11^. With the appropriate 
integrals in hand, the desired convergence essentially reduces to the fact that (1 + -)'^ — )■ e" 
as n ^ oo. So, before, we begin the proofs, we make some elementary observations about 
this convergence. If e~" < a; < 1, then < 1 + < 1, so: 

nlog 1 H < n = logx =^ - 1 H < 1 (3-11) 



n / n X \ n ^ 

If additionally it holds that (logx)^ < |, then: 

nlogfl + i^Unfi^ + of^^'lUlogx + of^^ 



\2- 



1 / logx\" ^ / (logx) , 

- 1 + -^ =1 + ^ ' (3.12) 
X \ n I \ n I 
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For the proofs of Lemmas 13.21 and 13.31 it is most convenient to use the normahzation 
Area(ii') = n and A = 1 (breaking from our previous convention). Thus n will be a positive 
integer, K will have area n, and we let Hk = ^k,i and Pk = PK,n- 

Proof of Lemma \3.^ Let X denote either N or A. In the derivation of the integral geometric 
expressions, we treat YIk and Pk simultaneously. 
The following formula is tautological: 

E[x] = / / ix{p,e)de (3.13) 

Jr/2tt JK 



dP{W{9)=p) (3.14) 

h=0 



Now let us derive expressions for Ix{p, 0) for the uniform and Poisson models respectively. It 
will be convenient to let y^^p equal /(p, 9)~^ times the distance from p to dK in the positive 
6 direction. For every angle 6', the coordinates Xp^ := exp(— 74(p, 9)) and Up ^ give a bijection 
between K and [e"", 1] x [0, 1]. It will prove very useful to express points in K in terms of 
these coordinates, mostly because dPiWuiO) = p) = exp{—A{p, 6)) dp = dx dy. 



(3.15) 
(3.16) 

(3.17) 
(3.18) 

And we also observe that dP{W{6) = p) equals: 

Forn^: e^^{-A{p,9))dp (3.19) 
(:-^)"_^., ,3.0) 

Using coordinates x and y in the integral f l3.13p and substituting our expressions for 
Ix{Pi d)-, we observe that we can integrate out yp^ in every case. The reader can check that 
the final expressions are: 

E[X] = f [ Ix{x,6)dxde (3.21) 



First, let us observe that 




h)\w{e)=p]l^^ 


is equal to 


For Uk 


and X = N: 


lylofip,or 




For Pk 


and X = N: 


\ylefiP.o?- 


n — 1 
- A{p, 9) 


For Uk 


and X = A: 


\ylef{p.o? 




For Pk 


and X = A: 


\ylef{p.of 
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where Ix{x,9) equals: 



For Uk and X = N: -f(p, df (3.22) 

6 

For Pk and X = N: \f{p, Of- (l + f 1 - -\ (3.23) 

6 X \ n J \ n J 

For and X = A: -f{p,ef (3.24) 



6* 



n-1 



For P;^ and X = A: -f(p, eV- 1 + (3.25) 

6 a; \ n J 

We now proceed to use the representations (13.211) and (I3.22p - p.25p to show that the expec- 
tations of X{IIk) and X{Pk) are the same up to a relative error of 0{n~^~^'^). 

First, observe that our estimate on the growth of / (Lemma 12.61) shows that cutting off 
the integral (13.211) to x > n^^ for some large fixed B incurs a relative error of no more 
than n~^~^'^. Now for x G [n^^ , 1], we may use (13.121) to see that the relative error incurred 
by replacing ^ (l + ^2££j" Iqj i jg j^q more than ^'"^"^ . Observe also that for x G [n~^ , 1], 
we know that replacing 1 + with 1 incurs a relative error of no more than These 
operations suffice to transform between the expressions for E[X(n/^)] and E[X(Pft:)], so we 
have shown that they are equal up to a relative error of 0(n^^+^). 

Thus to finish the proof, we just need to show that: 

E[X{U)]n~^+' = o(VVarX(n)) (3.26) 

By a result of Barany and Reitzner [2], VarX(n) ^ E[X(n)], so it suffices to show that 
A/E[X(n)] = o(n^"'). It is trivial to see that E[X(n)] < n, so we are done. □ 

Proof of Lemma lS.Sl This proof follows the same outline, so we will be a little less explicit; 
the interested reader can write down the long integrals if they so desire. Again, we let X 
denote either N or A. 

Our plan is to show that E[X^] is the same in the two cases up to a relative error of 

We think of X as being the integral of a random measure fix on M/27r. This random 
measure is just given by the family of variables X{a,P) (explicitly, the measure of the 
interval [a,/3] is X{a,/3)). Now X^ is just the total mass of fix ® /^x on (M/27r)^. Using 
linearity of expectation, we just need to take the d6 dip integral of the expectation of X{6,6 + 
d6)X{ilj,ip + dip). This expectation in turn, we condition on W{6) and W{ip), writing it as 
an integral dP{{W{6),W{ip)) = {p,q)) over K x K. We will often implicitly use the fact 
that the integrand is positive. 

The first step is to show that we may remove the region where either exp(— y4(p, 9)) < n^^ 
or exp(—74(g, '?/')) < and incur a relative error of 0{n~^^'^). By symmetry, let us deal 
with the region where exp(— A(p, ^)) < n^^ . Then the contribution to the total integral 
representing E[X^] is just: 



j I ±E[Xi9,9 + h).X\W{9)=p] 



/2tt {peK:A{p,e)>B\ogn} 



dP{W{9)=p)d9 (3.27) 

h=0 
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Now the X in the expectation contributes at most a muhiphcative factor of n + 2. With this 
X removed, the integral becomes something we aheady estimated in the proof of Lemma 
13.21 as being 0(n~^+'^). Thus we are done. 

The second step is to show that on the region where A{p, 6) and A[q, ip) are both < 
i^logn, the integrands (corresponding to E,[X{6,6 + d6)X{'ilj,tp + dip)] in the respective 
models) are equal up to a relative error of 0{n~^^'^). As before, this splits up into two 
problems: 

First, we need to show that the probability densities dP{{W{9),W{ip)) = {p,q)) in the 
cases of Uk and Pk are the same up to a relative error of 0{n~^^'^). As before, we may express 
dP{{W{9), W{i))) = {p, q)) elementarily in terms of A{p, q, 9, ip) := Area(C(p, 9) U C{q, ijj)). 
Then the fact that this quantity is O(logn) means we may apply (13.121) to see that the 
densities are equal up to a relative error of 0{n~^^'^) (note that this is true even for the 
singular part of the measure dP{(W{9), W{ip)) = {p, q)) occurring on the diagonal p = q). 

Second, we need to show that the incremental expectations E[X(6', 9 + d9)X{ip, ip + dip)] 
are the same up to a relative error of 0(n~^"^^). Again, this just involves writing equations 
such as (I3.15p -f l3.18p . Then the fact that A{p,q,9,ip) = 0(log?T,) shows easily that they 
coincide up to a relative error of 0{n~^^'^). Though it presents no difficulty in the proof, 
one should note that when X = N, there is a singular component to the measure E[X(6', 9 + 
d9)X{ip, ip + dip)] on the diagonal 9 = ip. 

We have shown that: 

\E[X{Uk)^] - E[X{Pk)^] I < n-^+'E[X{UK)^] (3.28) 

Now E[X{Uk)^] = E[X(nK)]^ + VaiXiUK). Thus we have: 

\E[X(nK)^]-E[X{PK)^] \ < n-i+^max(E[X(n^)]2,VarX(n;^)) (3.29) 

In the proof of Lemma \3.2\ we showed that E[X(ni^)] and E[X(Pi^)] are the same up to a 
relative error of 0{n~^^'^). This implies then that: 

\E[X{Uk)]^ - E[X{Pk)]^\ < n-^+'E[X{UK)]'' (3.30) 

Thus it follows that: 

iVarX(nA') - VarX(P^)| < max(E[X(ni^)]^ Var X(n/^)) (3.31) 

Thus to finish the proof, we just need to show that r?,^"'^^'^E[X(n)]^ = o(Var[X(n)]). By a 
resuh of Barany and Reitzner [2], VarX(n) > E[X(n)], so it suffices to show that E[X(n)] = 
o{n^~'^). It is a well known estimate (see |1]) that E[X(n)] ^ n^^^, so we are done. □ 
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